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SUMMARY 

An  experimental  and  theoretical  study  of  heat  transfer  and  flow  in  a 
rotating  duct  is  under  development  at  ONERA.  The  duct  simulates  in  a  simplified 
manner  an  internal  cooling  passage  of  the  blade  of  a  turbine  rotor. 

It  has  been  experimentally  demonstrated  that  the  rotational  velocity  causes 
an  increase  in  thermal  exchange  by  convection  in  the  duct.  This  phenomenon  is 
explained  by  secondary  flows  due  to  Coriolis  force. 

The  theoretical  aspect  of  the  work  involves  the  solution  of  the  Navier- 
Stokes  equations  in  three  dimensions  (time  marching,  ADI  method) .  For  both 
laminar  and  turbulent  flows,  the  results  obtained  demonstrate  the  significance  of 
secondary  structures  (having  both  two  and  four  vortices) . 

Also,  by  using  a  rotation-curvature  analogy,  visualisation  of  the  flow  with 
air  or  hydrogen  bubbles  in  a  water  duct  gives  rise  to  the  same  flow  structures. 
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1  INTRODUCTION 

Experimental  work  has  been  conducted  at  ONERA1,  with  the  aim  of  studying 
the  influence  of  the  effects  of  rotation  on  convective  heat  transfer  within  the 
internal  cavities  of  turbine  rotor  blades.  This  work,  in  common  with  that 
carried  out  with  other  experimental  installations2-4,  in  which  the  axis  of  the 
duct  under  consideration  is  perpendicular  to  the  axis  of  rotation,  illustrates 
the  differences  in  the  convective  exchanges  compared  with  those  obtained  in  a 
static  duct. 

These  differences  are  related  to  secondary  flows  caused  by  Coriolis  forces, 
also  to  the  effect  of  buoyancy  close  to  the  walls,  which  is  for  the  most  part  due 
to  the  centrifugal  force. 

In  parallel  with  this  experimental  approach,  there  has  been  a  theoretical 
one  in  hand  which  employs  the  numerical  modelling  of  flow  in  a  rotating  duct. 

This  approach  involves  the  solution  of  the  Navier-Stokes  equations  in  three 
dimensions.  A  number  of  programs  have  been  compiled.  The  first  of  these  can 
accommodate  only  simple  shapes  and  may  be  regarded  as  a  preliminary  stage  of 
development.  In  a  second  stage,  effort  has  been  directed  towards  a  more  general 
program,  which  is  being  validated.  The  interesting  feature  of  this  program, 
apart  from  predicting  the  effects  of  rotation  in  a  duct  of  simple  shape,  will  be 
the  representation  of  flow  and  heat  transfer  in  the  actual  internal  cavities  of 
rotor  blades  (Fig  l)5.  This  implies  the  simulation  of  effects  associated  with 
large  curvatures  (Fig  2)  or  disturbances  related  to  "turbulence  initiators" 

(rib3,  web3,  pin-fins.  Fig  3)  with  or  without  rotation. 

2  GENERAL 

As  pointed  out  in  the  Introduction,  the  internal  cavities  of  modern  turbine 
blades  are  of  complex  shape.  This  complexity  is  related  to  the  pursuit  of  an 
artificial  increase  in  the  turbulence  of  the  flow  by  means  of  shape  (successive 
curves,  turbulence  initiators)  to  lead  to  an  enhancement  in  convective  exchange. 

The  rotation  of  the  rotor  blades  adds  to  the  effect  of  the  shape. 

Although  the  combined  study  of  complex  shape  and  rotation  is  of  practical 
interest,  it  does  not  allow  the  separate  effects  to  be  readily  distinguished. 

The  normal  approach  is  to  study  each  of  the  phenomena  independently  before  com¬ 
bining  them. 

Furthermore,  the  prediction  of  these  turbulent  flows  (10000  £  ReDH 
<,  €0000),  which  can  be  carried  out  only  by  the  solution  of  the  Navier-Stokes 
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equations  in  three  dimensions,  requires  these  programs  to  be  validated  with  the 
model  of  turbulence  employed.  The  latter  approach  becomes  progressively  complex. 

As  a  first  step  therefore,  the  chosen  shapes  will  be  straight  ducts  with 
rectangular  sections. 

In  Fig  4,  a  simplified  internal  cooling  passage  of  a  turbine  blade  is 
shown,  together  with  the  Coriolis  and  centrifugal  forces  caused  by  rotation.  Two 
straight  "centrifugal"  and  "centripetal"  ducts  are  joined  by  a  curved  section. 

The  experimental  results  demonstrate  a  significant  change  in  the  transfer  of  con¬ 
vective  heat  with  rotational  velocity1-4.  An  overall  increase  in  convective  heat 
exchange  as  a  result  of  local  effects  is  observed;  an  increase  along  the  pressure 
faces  of  the  internal  duct,  and  a  reduction  along  the  suction  faces. 

These  changes  in  heat  transfer  with  rotational  velocity  are  explained  by 
the  effects  of  secondary  flow  associated  with  Coriolis  force  and  also  buoyancy, 
which  becomes  more  significant  with  temperature  difference  between  the  wall  and 
the  flow. 

The  secondary  flows  in  particular  are  responsible  for  a  distortion  of  the 
velocity  profiles,  compared  with  those  obtained  in  a  static  duct,  and  therefore 
also  of  the  differences  in  coefficients  of  friction.  The  resulting  local  vari¬ 
ations  in  the  exchange  coefficients,  associated  with  mixing  of  the  flow,  lead  to 
an  overall  increase  in  heat  transfer. 

The  equations  governing  the  flow  of  a  compressible  fluid  in  a  straight  duct 
are  given  below  for  a  system  of  axes  relative  to  the  duct  (Fig  5) .  They  include 
terms  relating  to  centrifugal  force  (pu>2x,  po>2y)  and  Coriolis  force  (2pti)v  and 
-2pti)u)  6 . 

We  have  equations  of: 


continuity 


H  *V-  '  <«> 


0  , 


motion 


+  V  .  (pv  V  +  pj 


J) 


pQ2x  +  2pv  A  Q  , 
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energy 


_a_ 

9t 


(pE)  +  V .  (pV  E  +  V  p  -  V  .  J  +  g  ) 


pQ2*  .  v 


ft  =  vector  rotation  (0,  0,  <o) 
n  =  radius  (x,  y,  0) 

V  =  velocity  vector  (relative)  (u,  v,  w) 

1  =  tensor  variable 

2  =  frictional  constraint  tensor 

g  =  heat  flux  vector 

In  the  dimensional  analysis  the  following  non-dimensional  quantities  are 


significant : 

Nusselt 

Nu 

=  hDH/k 

characterising 

the 

heat  transfer 

Reynolds 

Re 

=  pGDH/p 

characterising 

the 

flow  turbulence 

Rossby 

Ro 

=  <odh/u 

characterising 

the 

Coriolis  aspect 

Ekmann 

Ek 

2  , 

=  Re  .  Ro  *  pwD  /p 

H 

characterising 

the 

Coriolis  aspect 

Rayleigh 

Ra 

2  2  3 

=  V  (0  HD  (Tp  -  Ta  .  ) 

H  1 

.  Pr  characterising 

the 

centrifugal  aspect 

These  quantities  are  employed  in  order  to  group  together  the  results 
relating  to  the  influence  of  rotation  on  the  heat  transfer. 

3  THERMAL  MEASUREMENTS 

A  general  view  of  the  test  rig  is  given  in  Fig  6.  Its  dimensions  and 
characteristics,  as  well  as  the  experimental  conditions,  were  chosen  to 
facilitate  the  taking  of  measurements.  Also,  in  order  to  use  as  simple  a 
technology  as  possible,  the  rotational  velocity  and  the  temperature  of  the  test 
models  were  reduced  by  comparison  with  the  values  met  in  an  actual  blade 
(t0max  =  5000  rev/min  and  Tpmax  -  400°C  ,  instead  of  13000  rev/min  and  1000°C 
for  a  rotor  blade) . 

Two  similar  test  models  (Figs  7  and  8)  are  able  to  turn  inside  an  electric 
heater  (Tf-ouR.max  “  1200°C  and  P  -  45  kW)  .  This  assembly  is  placed  inside  a 
vacuum  chamber  (pressure  <20  Pa) .  Thus  the  models  are  externally  heated  by 
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radiation  from  the  electric  heater;  external  convection  is  negligible  on  account 
of  the  low  pressure  in  the  chamber.  They  are  cooled  by  the  circulation  of  air 
entering  at  the  bottom  of  the  rotating  shaft  and  leaving  at  the  top;  the  flow 
rate,  in  the  range  5-50  g/s,  is  measured  by  an  ultrasonic  flowmeter. 

Four  sets  of  rotating  seals  ensure  the  absence  of  leaks  in  the  air  supply 
and  the  vacuum  chamber. 

An  electric  motor  drives  the  shaft  supporting  the  models  via  a  toothed 

belt. 

Measurements  are  made  using  thermocouples  attached  to  the  various  sections 
(Fig  7) .  An  infrared  pyrometer  of  short  time-constant  (At  -  10  pis)  and  spatial 
resolution  of  the  order  of  1.2  mm  allows  contours  of  the  radiation  temperature  to 
be  obtained  in  a  section  of  the  test  model7.  In  Fig  8,  an  example  of  a  contour 
obtained  by  pyrometer  is  given,  with  a  model  similar  to  the  cooling  passage  shown 
in  Fig  4.  From  two  angles  of  view  (8  -  ±4  5°)  the  temperature  contour  is  con¬ 
structed  for  two  rotational  velocities.  Distortion  of  the  contour  with  increase 
in  rotational  velocity  may  be  clearly  seen.  The  points  c,  d,  e  and  f  define 
the  boundaries  of  a  section  (x  -  r$) /DH  -  9.5  centrifugal  duct,  and  it  is  to  be 
noticed  that  the  transfer  of  heat  by  convection  is  greatly  increased  for  the 
pressure  face  (the  concave  surface  of  the  rotor  blade) ,  since  the  wall  tempera¬ 
ture  is  less,  compared  with  the  case  of  low  rotational  velocity  (0)  -  50  rad/s) . 

On  the  other  hand,  for  the  suction  face  (the  convex  surface  of  the  blade)  the 
exchange  is  less.  The  significance  of  the  measurements  carried  out  by  infrared 
pyrometry  is  the  direct  observation  of  the  thermal  effects  asociated  with 
rotation  and,  in  particular  of  the  Coriolis  effect,  which  gives  rise  to  flow 
phenomena  which  are  different  at  the  opposing  concave  and  convex  faces.  The 
differences  noted  between  the  values  obtained  from  the  thermocouples  may  be 
explained  by  reflections  of  the  radiation  from  the  heater  on  to  the  part  of  the 
surface  being  observed.  The  output  of  the  pyrometer  is  of  the  form: 


u 


'VW'-V 


L,  )  dl  . 

A,env 


internal  function 
of  the  pyrometer 


actual  emission 
from  the  surface 


reflection  from 
the  surroundings 
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The  exchange  coefficients  are  derived  by  carrying  out  a  thermal  balance  at 
each  thermocouple  under  both  transient  and  steady-state  conditions: 


pee 


dTp 

dt 


E<ginc  -  "  h<TP  "  T°>  +  9cond  • 


For  the  centrifugal  duct  of  the  model  shown  in  Fig  7,  the  ratio  of  the 
Nusselt  number  with  rotation  (Nuq)  to  the  Nusselt  number  in  a  static  condition 
(Nu0)  is  given  as  a  function  of  the  Rossby  number  which  describes  the  Coriolis 
aspect.  The  values  correlate  quite  well  with  the  Rossby  number,  and  as  stated 
above,  the  convective  exchange  is  increased  significantly  at  the  concave  face  and 
in  this  case  is  decreased  to  a  smaller  extent  at  the  convex  face.  This  result  is 
consistent  with  that  obtained  by  Morris  2  and  Johnson3. 

4  NUMERICAL  SOLUTION  OF  THE  NAVIER-STOKES  EQUATIONS  FOR  A  ROTATING 

DUCT 

The  prediction  of  flows  and  heat  transfers  in  an  internal  cooling  passage 
of  a  rotor  blade  requires  the  solution  of  the  Navier-Stokes  equations  in  three 
dimensions.  This  modelling  mu3t  be  validated  for  a  certain  number  of  elementary 
configurations  before  being  used  for  any  shape.  The  comparison  of  computation 
with  the  experimental  results  given  in  section  3  is  an  example. 

Some  different  programs  were  developed  to  solve  the  equations  given  by 
system  1  of  section  2.  The  first  ones  allowed  flow  in  simple  shapes  to  be  dealt 
with,  in  this  case  a  rotating  straight  duct  or  loop.  The  principle  behind  the 
computation  is  normal,  and  comprises  a  centred  finite  difference  method,  based 
upon  a  three-dimensional  curvilinear  orthogonal  system,  whose  conservative  nature 
is  assured  by  the  integration  of  the  continuity  equations  over  the  elementary 
volumes  associated  with  each  node  of  the  network.  These  volumes  of  integration 
are  formed  inside  the  domain  of  computation  from  three  consecutive  points  in  each 
direction.  In  the  limits,  these  imaginary  "sliding  velocities"  are  stored,  and 
the  continuity  equations  are  integrated  over  the  half-volume  associated  with  the 
point  under  consideration,  taking  account  of  the  limiting  conditions  through  the 
medium  of  the  wall  functions.  The  spatial  discretisation  provides  a  system  of 

5  Nx  Ny  Nz  equations  of  evolution  (Nx>y  and  z  -  number  of  points  in  the  network  in 
the  different  directions),  which  is  integrated  in  the  course  of  time  by  a  semi- 
explicit  method.  The  inversion  operators  are  linearised,  then  factorised  in  each 
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direction  (by  the  ADI  method) ,  and  the  computation  consists  of  a  series  of 
inversions  of  tridiagonal  matrices  by  (5  x  5)  blocks8. 

At  the  moment  effort  is  directed  towards  a  more  general  program  concerned 
with  both  the  shape  (of  any  curvilinear  network)  and  the  limiting  conditions. 

The  technique  of  integration  in  time  is  analogous  to  that  just  described.  The 
spatial  discretisation  employs  a  method  of  finite  volumes  with  flux  splitting  (of 
the  Steger  and  Warming  type9)  and  a  diagram  decentred  at  five  points. 

The  principal  feature  of  this  program,  which  was  conceived  for  the  predic¬ 
tion  of  flow  in  a  diesel  engine,  is  the  ability  to  modify  the  shape  at  each  3tep 
in  time.  With  regard  to  the  mesh,  a  transformation  which  changes  the  physical 
domain  of  computation  into  an  assembly  of  patalellepipeds  facilitates  the  appli¬ 
cation  of  the  method.  The  limiting  conditions  are  very  diverse:  they  include 
temperature  or  flow  at  areas  of  the  wall,  and  local  injection  or  suction  of  the 
flow.  For  the  accuracy  achievable  with  this  program  the  reader  should  refer  to 
Refs  10  and  11. 

In  Fig  9  an  example  of  the  results  of  asymptotic  three-dimensional  compu¬ 
tation  in  laminar  flow  is  given.  The  flow  conditions  are:  ReDH  *=  1000  ; 

Ro  «  0.125  and  EK  *  125  .  and  this  computation  corresponds  to  a  situation  in 
which  x  is  much  greater  than  a  and  b  .  The  longitudinal  velocity  gradients 
du/dx  are  assumed  negligible.  The  quantity  (1/p)  (dp/dx)  -  c^x  is  fixed;  it  is 
homogeneous  over  the  entire  section,  and  it  determines  the  rate  of  flow.  Fig  9a, 
which  shows  a  system  of  two  vortices,  corresponds  to  a  40  x  40  mesh  and  was 
obtained  after  1000  iterations  (Dh  “  10  mm,  t  *  5.10_5s,  CFL  *  80).  Fig  9  shows 
another  case  with  Reno  *  500  ,  Ro  »  0.25  ,  and  EK  »  125  ,  the  secondary  flow 
structure  being  changed  into  a  system  of  four  vortices.  The  intense  recircu¬ 
lation  of  the  two  supplementary  vortices  is  to  be  noted.  This  type  of  result  is 
in  agreement  with  that  obtained  by  Khesgi  and  Scriven1?,  or  Speziale  et  al13. 

An  example  of  the  results  obtained  in  three  dimensions  in  a  rotating  loop 
is  given  in  Fig  10.  The  (16  x  32)  x  96  mesh  represents  a  duct  of  rectangular 
section  (a/b  “  2)  comparable  to  the  one  investigated  and  shown  in  Fig  7.  The 
recirculations  and  the  axial  velocity  contours  are  illustrated  in  Fig  10.  They 
demonstrate  the  complexity  of  the  flow  which  in  this  instance  is  turbulent: 

Re  ”  40000,  Ro  ”  0.06  and  Ex  -  2400  .  The  inlet  conditions  are  a  specified 
stagnation  pressure  and  enthalpy.  Downstream,  only  the  static  pressure  is 
specified.  Flow  is  zero  at  the  walls  (adiabatic  flow) . 
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The  results  obtained  with  the  simplified  computer  program  (first  stage,  see 
above)  demonstrate  the  feasibility  of  the  confutation  for  turbulent  conditions. 
They  are  encouraging,  and  the  numerical  study  will  be  followed  by  the  program 
under  development  (second  stage),  which  will  then  allow  a  quantitative  analysis 
of  the  thermal  exchanges,  using  different  walls. 

5  VISUALISATION  OF  FLOW  LINKS  IN  A  CURVKD  DUCT 

An  experimental  approach  to  the  visualisation  of  secondary  flows  in  a 
curved  duct  was  embarked  upon.  In  fact,  analysis  of  the  equations  of  motion 
demonstrates  the  similarity  of  behaviour,  with  respect  to  centrifugal  effects,  of 
flow  in  a  straight  rotating  duct  and  that  encountered  in  a  curved  duct.  In  both 
cases,  the  trajectories  of  the  particles  in  a  Galilean  frame  of  reference  are 
curved  and  the  phenomena  are  determined  by  interaction  between  the  inertial 
forces  due  to  the  curvature  and  the  non-uniformity  in  flow  velocity  due  to  wall 
friction14 . 

The  significance  of  such  visualisation  is  to  give  prominence  to  the  second¬ 
ary  structures  brought  into  being  by  this  effect,  also  to  throw  light  on  the  two- 
and  four-vortex  duality,  to  obtain  information  on  the  thickness  of  the  recircu¬ 
lation  zones,  and  hence  that  of  the  vortex  centres  themselves. 

In  practical  terms,  the  rotation-curvature  analogy  is  achieved  by  equating 
the  Coriolis  and  centrifugal  accelerations  associated  with  the  curvature: 

2 

20)u  =  u  /R  =*  R  •  u  /2( 0  , 


and  in  terms  of  the  equivalent  Rossby  number  for  the  curvature: 


Ro 


courbure 


Dh/2R 


The  air-water  analogy  is  achieved  by  conserving  the  Dean  number 


Dn  -  Reno  .  V  1  with  y  =  2R/D  -  1/Ro 

DH  T  T  H  courbure 
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The  values  chosen  for  the  experiment  described  here  are  (Fig  11): 


R  “  0 . 12  m  ;  D|]  ■  40  nun  ;  Re^^  ■  18000  ;  ^^courbure  "  0.166  , 


with  a  duct  of  rectangular  section  having  transparent  Altuglass  walls,  of  height 
a  -  60  mm  and  width  b  -  30  mm  (a/b  «  2,  Fig  7) .  These  dimensions  were  chosen 
so  as  to  allow  practical  visualisation  of  the  flow  using  bubbles  of  air  or  hydro¬ 
gen.  Hydrogen  bubbles  are  obtained  from  two  grids  made  up  of  0.5  mm  diameter 
wires  in  a  2  x  2  mm  mesh;  they  are  spaced  2  mm  apart,  and  have  40  V  applied 
across  them.  In  a  vertical  position  they  occupy  almost  the  whole  cross-section 
of  the  duct.  Injection  takes  place  12.5  x  Djj  (500  mm)  from  the  curved  part. 

Illumination  of  a  section  2  mm  thick  allows  the  flow  lines,  and  therefore 
the  reconstruction  of  the  secondary  flow  (Fig  12),  to  be  visualised  with  equal 
ease  both  vertically  and  horizontally.  In  the  first  case  the  illumination  is 
arranged  to  follow  the  curvature  of  the  duct;  in  the  second,  it  is  plane. 

In  Fig  12.1,  the  flow  is  visualised  in  the  region  of  the  upper  wall  and  the 
circulation  is  centripetal;  the  recirculation  then  follows  the  inner  wall 
(Fig  12.2).  In  Fig  12.3  (with  horizontal  illumination  in  the  middle  of  the  duct) 
a  slight  centrifugal  effect  appears.  Along  the  external  wall  (Fig  12.4.1)  the 
two  centres  stand  out  clearly  with  an  increase  in  the  flow  in  the  upper  half  and 
a  fall  in  the  other  (see  also  the  recirculations  in  Fig  10) .  An  instability  at 
the  external  wall  in  the  central  area  is  evident  on  examining  Figs  12.4.1  and 
12.4.2.  In  this  separation  zone  of  the  two  vortices,  the  flow  rises  or  falls 
intermittently;  this  may  be  explained  by  the  appearance  and  disappearance  of  two 
small  supplementary  vortices.  In  Figs  12.1  to  12.4.2  visualisation  was  achieved 
by  air  bubbles  and  the  flow  may  be  considered  established  at  the  entrance  of  the 
curved  duct,  as  the  straight  part  of  rectangular  cross-section  is  1  m  long,  ie 
25  x  Dh . 

Injection  of  hydrogen  bubbles  was  employed  in  order  to  obtain  bubbles  of 
smaller  diameter  (d  s  5/100  mm)  than  those  of  air  (d  o;  0.5  mm),  and  which  there¬ 
fore  follow  the  flow  better.  There  is  in  fact  no  noticeable  difference  in 
behaviour  of  the  two  types  of  bubbles  for  flow  velocities  greater  than  0.4  m/s. 

However,  by  inclining  the  system  for  generating  the  hydrogen  bubbles  by  a 
few  degrees,  a  structure  having  four  vortices  became  evident  (Fig  12.4.3,  flow 
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visualisation  near  the  outside  wall) .  With  a  significant  angle  of  inclination 
(45°  in  Fig  12.4.4)  the  structure  becomes  asymmetrical. 

Thus,  the  natural  structure  appears  to  be  a  two-vortex  system;  however,  a 
small  disturbance  can  induce  a  four-vortex  system.  A  more  significant  disturb¬ 
ance  leads  to  asymmetric  circulations.  These  results  are  not  changed  qualitat¬ 
ively  when  the  Reynolds  number  varies  in  the  range  10000-60000. 

In  Fig  12.1  it  is  to  be  noted  that  the  path  of  the  bubbles  makes  an  angle 
of  about  45°  with  respect  to  a  tangent  to  the  axis  of  the  duct.  For  the  internal 
and  external  faces  (Fig  12.2  and  12.4.1)  the  angle  between  the  horizontal  and  the 
path  of  the  bubbles  is  smaller,  25°  and  15°  respectively. 

These  angles  become  smaller  as  the  Reynolds  number  increases. 

The  zones  of  recirculation  are  thin  and  of  the  order  of  5-10%  of  the 
hydraulic  diameter. 

It  follows  that  correct  modelling  under  turbulent  conditions  in  this  type 
of  flow  may  De  carried  out  only  by  refining  the  mesh  adjacent  to  the  walls. 

In  conclusion  to  this  section,  a  question  of  a  purely  practical  nature 
arises:  what  is  the  effect,  where  convective  exchanges  are  concerned,  of  ‘-hese 
different  structures,  ie  two-vortex,  four-vortex,  or  asymmetrical? 

In  fact,  in  addition  to  the  phenomena  described  above,  entry  effects 
playing  a  significant  part  in  the  development  of  the  secondary  structures  in 
actual  situations  (such  as  the  internal  cooling  passage  of  a  blade)  may  appear 
and  render  any  predictions  difficult. 

6  CONCLUSION 

Knowledge  of  the  flow  and  heat  transfer  in  the  internal  ducts  of  the  rotor 
blades  of  turbines  is  necessary  in  order  to  optimise  the  shape  of  these  cavities 
and  obtain  the  best  thermal  efficiency. 

At  ONERA  a  test  rig  has  been  set  up  to  evaluate  the  influence  of  the 
effects  of  rotation  on  the  convective  exchanges  in  a  duct  which  simulated  in  a 
simplified  manner  an  internal  cooling  passage  of  the  rotor  blade  of  a  turbine. 
Experiments  have  been  carried  out  using  two  models  comprised  of  tubes  having 
rectangular  sections  of  ratio  2  and  0.5  respectively  between  the  sides.  Using 
these  shapes  the  coefficients  of  convective  exchange  were  significantly  modified 
when  compared  with  those  encountered  using  a  static  duct.  When  the  rotational 
velocity  was  increased,  the  overall  convective  exchange  increased;  there  was  a 
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significant  increase  on  the  pressure  face  of  the  duct,  and  a  decrease  along  the 
suction  face. 

These  effects  may  be  largely  explained  by  secondary  flows  induced  by 
Coriolis  force. 

In  parallel  with  this  approach,  numerical  modelling  by  solution  of  the 
Navier-Stokes  equations  in  three  dimensions  is  under  development.  The  results 
obtained  for  laminar  flow  are  in  accordance  with  those  published  by  other 
authors12'13,  and  for  certain  values  of  Rossby  and  Ekmann  numbers  recirculations 
appear  with  either  two  or  four  vortices. 

For  turbulent  flow,  a  complete  three-dimensional  modelling  of  a  loop  with 
centrifugal,  curved  and  centripetal  ducts  has  been  carried  out.  The  complexity 
of  the  secondary  flows  is  very  apparent,  especially  in  the  curved  part;  distor¬ 
tion  of  the  axial  velocity  contours  is  evident,  as  predicted,  to  a  considerably 
less  degree  than  in  the  laminar  mode.  This  kind  of  computation  will  be  continued 
using  a  Navier-Stokes  three-dimensional  program  undergoing  validation,  which  will 
allow  quantitative  analysis  of  the  thermal  exchanges  with  the  walls. 

Secondary  flows  which  arise  in  a  straight  rotating  duct  may  be  studied 
experimentally  and  more  simple  in  a  curved  duct  by  effecting  a  rotation-curvature 
analogy.  The  experimental  conditions  are  defined  from  a  situation  which  has  been 
analysed  on  a  thermal  basis  (model  a/b  =  2,  ReDH  =  18000  and  Ro  -  0.16). 

The  visualisation  of  the  paths  of  air  or  hydrogen  bubbles  in  the  curved 
duct,  using  illumination  in  different  horizontal  and  vertical  planes,  demon¬ 
strates  a  system  of  two  vortices.  However,  a  moderate  flow  disturbance  intro¬ 
duced  at  12.5  x  dh  upstream  of  the  curved  duct  causes  this  structure  to  be 
changed  into  a  four-vortex  system.  A  greater  disturbance  gives  rise  to  an 
asymmetric  situation. 

These  qualitative  results  agree  with  the  predictions  of  Soh16  for  laminar 
flow,  in  respect  of  the  influence  of  upstream  conditions  on  secondary  structures. 
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C,  Cp  specific  heat 

Dh  hydraulic  diameter 

e  wall  thickness 

E  energy  -  CVT  +  V2/2 

Ef(  Ekmann  number 

Fce  centrifugal  force 

Fco  Coriolis  force 

h  heat  exchange  coefficient 

H  eccentricity  =  (re  +  L/2> 

k  thermal  conductivity 

Lx  monochromatic  radiation 

Lx°  black  body  radiation 

m  mass  flow  of  air 

Nu  Nusselt  number 

P  pressure 

Pr  Prandtl  number 

Prt  turbulent  Prandtl  number 

qCond  conduction  heat  flux 

qinc  incident  radiation  heat  flux  from  heater 

r,  R  radius 

Ra  Rayleigh  number 

Re  Reynolds  number 

Ro  Rossby  number 

T  mean  temperature 

Ta  air  temperature 

TFour  heater  temperature 

Tp  wall  temperature 


t 


time 
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u,  v,  w  components  of  flow  velocity 

.  'A 

V  velocity  =»  (u  +  v  +  w) 

x,  y,  2  co-ordinates 

Greek  symbols 

P  coefficient  of  expansion 

e  emissivity 

\  wave] sngth 

(1  dynamic  viscosity 

V)  kinematic  viscosity 

p  density 

a  Stefan's  constant 

<o  rotational  speed 

Suffixes 

o  static  duct 

e  entrance  of  straight  duct 

i  inlet  conditions 
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DISCUSSION 

D.K.  Hennecke,  Technische  Hochschule,  Darmstadt,  Germany: 

In  your  measurements,  what  kind  of  velocity  distribution  in  the  inlet  cross 
section  of  the  test  section  did  you  try  to  achieve  and  how  did  you  do  it? 

Author's  reply: 

The  shape  of  the  inlet  part,  just  before  the  test  section  (centrifugal 
duct)  has  been  chosen  with  the  industrial  partner  (SNECMA)  in  order  to 
have  a  similar  shape  of  the  channel  before  the  blade. 

Therefore,  in  order  to  compute  accurately  the  flow  in  3D  (and  consequently 
the  friction  and  the  heat  transfer  coefficients  within  the  test  section)  it 
is  necessary  to  start  the  computation  near  the  axis  of  rotation  owing  to 
the  Coriolis  force  which  induces  secondary  flow  from  the  hub. 

J.  Moore,  Virginia  Polytechnic  Institute,  USA 

The  orientation  of  your  bend  is  at  90°  to  that  used  in  our  paper.  Are  you 
also  planning  to  test  different  bend  orientations? 

Author's  reply: 

With  the  first  test  section  <a/b  «  2),  the  objective  was  specifically  the 
study  of  the  two  channels: 

-  the  centrifugal  channel  one 

-  the  centripetal  channel  one. 

Therefore,  the  radius  of  the  curved  part  which  connects  the  two  ducts  has 
been  chosen  in  order  to  minimise  the  perturbation  of  the  flow  induced  by 
the  bend. 

But,  with  the  second  test  section  (a/b  -  0.5)  the  geometry  is  close  to  a 
real  cavity  of  blade  and  the  orientation  of  the  bend  is  at  90°  by 
comparison  to  the  first  test  section  and  like  the  one  computed  in  the 
paper  of  Mr  and  Mrs  Moore. 
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Fig  1  Example  of  the  complex  internal 
cavity  of  a  turbine  rotor  blade 


Fig  2  Example  of  a  simple  reversal  -  visualisation 
by  means  of  a  solution  of  Teepol  and 
water  -  ReDH  -  15000 


Figs  3,4&5 


Fig  3  Example  of  a  disturbance  -  duct  of  square  section  - 
visualisation  by  air  bubbles  <  RepH  - 15000 


Fig  4  Example  of  Internal  cavity 


Fig  5  Representation  of  a  straight  duct 
In  the  rotating  axis  system 
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Figs  6*7 
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Analogies:  rotation  =>  curvature 
air  =>  water 


Re  -  1800C;  Ro»  0  16 


Fig  12  Visualisation  of  lines  of  flow  In  a  curved  duct 


